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Abstract 

The condition for potential description of the wake waves, generated by flat or cylin- 
drical driving electron bunch in cold plasma is derived. 

The two-dimensional nonlinear equation for potential valid for small values of that 
is obtained and solved by the separation of variables. Solutions in the form of cnoidal 
waves,existing behind the moving bunch at small values of vertical coordinate, are ob- 
tained. In particular, at some boundary conditions, corresponding to blow-out regime in 
the underdense plasma,the solution represents by a solitary nonlinear wave. 

Approximate solution is also obtained using the method of multiple scales. 

The indications are obtained that the dependense of the amplitudes on longitudi- 
nal coordinate determines essentiallyeven in the first approximation,by driving bunch 
charge distribution. The wake wave amplitude can increase at some conditions along the 
longitudinal distanse from the rear part of the bunch. 

1 INTRODUCTION 

Analysis of the one dimensional longitudinal, transverse and coupled transverse lon- 
gitudinal plain nonlinear waves in cold relativistic plasma are given in the review [|l j 
(see therein the references on original works). One dimensional nonlinear longitudinal 
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waves,generated by the driving bunches with the infinite transverse dimensions, were con- 
sidered in 

In the present work the two-dimensional nonlinear wake waves,generated by the flat 
or cylindrical electron bunch, are discussed. 

The corresponding linear problem was considered in []9|, ||10|| and was found, in cor- 
respondance with previous result |Tl|,that the magnetic field in wake wave in linear 
approximation is equal zero. This result connected with the absence of the energy flow 
in the wake wave and the absence of the vortexes in plasma electron motion in linear 
approximation. 

In the one dimensional nonlinear treatment P]-[§] the magnetic field in the wake wave 
is also zero (by construction), due to the symmetry of the problem relative to transverse 
displasments. 

In two-dimensional wake wave the magnetic field is zero only when vortexes connected 
with the plasma electron motion are zero,which in this case is an additional requirement 
on the type of the motion of the plasma electrons. The wake waves in this case are po- 
tential, i.e. electric field components E z ,E y can be expressed as a gradient of one scalar 
function ip(y, z),and components of the plasma electrons current also can be expressed 
through one scalar function ip(y, z) 

The approximate nonlinear equation for potential (p can be obtained, using Maxwell 
equations and approximate equations of the motion. Equation for potential has an exact 
solutions with the separated variables for small values of the transverse coordinate. 

Among the solutions, which are finite, nonlinear waves by cnoidal nature, at some 
boundary conditions, associated with the blow-out regime, there exists the solution (on 
separatirix) in form of the solitary wave. 

For arbitrary values of the transverse coordinate the approximate solution for poten- 
tial ,using multiple scales perturbative method is found. 
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2 VORTEX-FREE WAKE WAVE 



Consider the wake wave generated in the cold neutral plasma,with the immobile 
ions,by the flat electron bunch, which has horisontal dimensions 2a much larger than 
vertical dimension 26,longitudinal dimension is 2<i.The charge density in the bunch is 
rib, electron plasma density is no, and we consider both overdense and underdense regimes. 

Bunch is moving along z-axis with the constant velocity v o < c in lab system. All the 
physical quantities in the question are considered as a function of vertical coordinate y 
and z = z — t> t.An electrical field,generated by the bunch \E X \ <C E y ^ 0, E z ^ and 
magnetic field B z = 0, \B y \ < \B X \ = \B\ ^ 0. 

Introduce the dimensionless variables and arguments by 



E = JlimmvlE' = ^*E> (1) 



B= \IAimrm%B>=^^B> 



z,y=kz,ky,k = — = -,n 6 =— ,n = — (2) 

Vq ttiVq n n 

where n is the arbitrary electron density, which is convinient to choose equal = in 
the underdense (n& > n ) case and n = no in the overdense (n& < n ) case. 
Following 113^1, P| introduce BFTCh-transformation of the variables 

< = ^ = ^ + V.); (4) 

(N -> n , /3 e2 — ► 0, when z, j/ -> +oo; v ex = 0). 

The Maxwell equations then can be rewritten in the following form (superscript prime 
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is ommited in what follows): 

(a) 



8(B + 0E S ) = 
dz 



(c) 



(5) 

d((3B + E y ) _ dE^ 
dz dy 



(d) ^ + ^ = (n -n b )-N(l + V z ) 

The continuity equation = ^ follows from (§.d),(§.a),(§.b). Using (§.a), (§b),(|.c) 
we have 

^ + (1 - /5 2 )|^ = rot x (J3NV) + rot x 0n b ) (6) 

which means that the magnetic field is zero in plasma (linear or nonlinear) wake wave 
only when 

rot(f3NV) = (7) 

i.e. the plasma electrons motion is vortex-free. 

In the following,we consider the region of the space, occupied by wake wave i.e. z < —d 
.Maxwell equations @ for wake waves under condition (|7D can be obtained putting in 
(§5 = and n b = 0. 

(a) ?*—NV m 



(8) 



(c) 



dE y _ dE z 



dz dy 



(d) ^ = n -N 
dy 



Then from (Re) follows that 



E = —gradf (9) 



i.e. the wake fields under condition (|7p,as it must be, are potential. 



3 THE BASIC EQUATION FOR THE POTENTIAL. EXACT SOLUTION IN SEPARA- 
BLE ARGUMENTS 

Consider Maxwell equations (||) for wake waves, when z < —d. From (^.a),(§[d) and 
© we have 

From hydrodynamic equation of the plasma wake wave electrons motion,using (|1]),(|2]),(|3|) 
it is possible to obtain the relativistic equation of motion for the V z component of the 
generalized velocity: 



-dV z , . dV, 



+ v v —^ = -w 1/2 



dz dy 



E z ^l + 2V z + ^j +(3 2 V z V y E y 

W = 1 + 2V Z + ^- {3V 2 ; 



(12) 



Neglecting terms with the squares of generalizid velocitycompared to the terms with 
the first power of that, the expression ( |12"D converted to 

dV z 



dz 

The solution of this equation, using ©,is 



E Z (\ + 3V S ) (13) 



V, = l(e- 3 * - 1) « -tp (14) 

with the condition ip = 0,when V z = 0. 

Substituting (|T4]) in ([0]) we have the basic equation for <p 

d 2 (p d 2 p> 

^ + ^ + ^ = (15) 

Nonlinear term in eq. ([15]) is proportional to |^ and can be large. In the linear approx- 
imation solution of eq. (fL5[) lost the y-dependence and describes the harmonic oscilation 



with the plasma frequency on z ^-dependence of the solution comes from boundary con- 
dition at z = —d and coincides with it for all z < — dThis is always the case,when wake 
waves are described as the product of the two functions from separate arguments y and 
z. Such a situation takes place in linear approximation [|3|j , p"T|j , |0| . 



The eq. ([15]) permits to search the solution with the separable variables 



<p(v,z)=<Pi(v)<P2(z) (16) 
if 2 + n <f2 



, - -<£ = -k (17) 
where k is a separation constant. The equations for ipi and ip 2 are: 

<fi" = k (18) 

ip 2 " + n ip 2 + kip 2 2 = (19) 



Due to the symmetry of the problem the solution of equation (|T8| ) must be symmetric 
on y; (fi(y = 0) = due to E y — at y — O.The solution of the linear problem @,|10| is 



concentrated in the region of the " trace" , falling outside it exponentially. Adopting the 
same picture of the potential flow for considering case too,the solution of eq. flT8] ) is 

which is valid for small values of y.It means that, the solution of eq. (^) in separable 
arguments exists only for small values of y < b; 

The equation ( |19|) is the equation for nonlinear oscilator,with nonlinear part of the 
force proportional to ip\ (for mathematical pendulum the first nonlinear term is propor- 
tional to ip\ see e.g. |TB| ). The general solution of this equation is given in the implicit 
form by 

-(*+rf)=±r — *» m . (20) 

sign ± corresponds to positive or negative subsequently and h is an energy con- 
st ant, defined by 

h = \p'\ + + (21) 
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and determined from boundary condition at z 
The function F(ip 2 ) is 



-d,<p 2 (-d) = <p ,<p' 2 (-d) = <p' 



(22) 



k 

2 rz ' 3' 

The separatrix,/i s = gp^o, is the tangent to F{ip 2 ) at its maximum pomt;F((p 2 ) has three 
real roots: double root equal zero and one root at B = ^ The roots of the equation 



h - FM = 0, 



(23) 



are q(i = 1,2,3) The different solutions (p0|) of the equation (|19D defined by the value 
of the /i, which in turn, depends on the boundary values ipo and <p' Q . Finite solutions for 
k < is C\ < ip 2 < c 2 , which is existed,when h < h s , i.e. c s < C\ < 0, where 



c, = - 



2\k\ 



and < C2 < c„ 



no 



(c m corresponds to the local maximum of the function F((p 2 ), equal to h s ).The third 
root of the eq fl23]) is c 3 : c m < c 3 < B 



In^.For A; > finite solutions exist, when 

2|fc| ' 



c 2 < ^ 2 < c 3 , < /i < /i s ; for k = 1, < Cl < _ nQ) -n < c 2 < 0, < c 3 < ^; 

For the cylindrical bunch with length 2c? and radius Ro eq. (|9]) is written in the form 



i a / % N 

r or \ or . 



V z 



The approximate equation of motion for V z has the same form as (|12"D and \4 approx- 
imately is equal to ip. Basic equation for <p is then 



#V 1 9 ( d<p\ 

dz^ + LP 2d-r VTr +n0Lf 







Solution of this equation in separable arguments <p(r, z) 
by solving the equations: 



tpi(r)ip 2 (z) can be obtained 



id_ (dcpA = k 

2 dr \ dr I 



d 2 (f 2 
dz 2 



+ n ip 2 + kcpl = 
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The last equations for (p 2 (z) coincides with the eq (|TJ) for flat case. Equation for (pi(r) 
has the solution finite for small r 



. Hence the cylindrical bunch case is described practically by the same equations as a 
flat one with the evident changes from (pi(y) to <pi(r). 

4 BOUNDARY CONDITIONS 

The definition <p = tpi(y)ip 2 (z) permits the transformation ipi = \k\ifi,ip 2 = | A;| _1 (^2-Then 
eqs. ([18]) and (plj) for (p\ and p> 2 will have the following forms: 



#' = ±l,#[ + n £ 2 ±$ = 0, (24) 

which corresponds to value of k = ±l.It means that separation constant k is arbitrary 
and can be choisen as k = ±l.In what follows the sign "bar" over (p\ j2 is ommited and k 
is choosen plus one,which provide meaningful results for our case, when z < —d; k = — 1 
is suitable for z > O.Then 

<p = <Pi(y)<P2(z) = {a + y 2 /2)<p 2 (z) (25) 

If at z = —d the physical quantities are 

E d z = E z (y, z = -d) = E z0 + E z2 y 2 = - (^j , (26) 
V z d ee V z (y, z = -d) = V z0 + V z2 y 2 = -<p(y, z = -d), 



the unknown constants (po, f'o, a> entering in the solution (TH 



arc 



From (|.a) and (ITJ 



<p' = -2E z2 , <p = -2V z2 , a = ^ = ^- (27) 

&z2 £ V z2 



Q 2 

NV Z = -^j = (pi(y)<p" 2 (z) = <Pi(y)(-n (p 2 - (pi) (2? 



and, when z=-d 



{NV z f = -(a + y 2 /2)(p + n )^ = V d & + n ), 



i.e. 



N d = ^ + n Q ,N d = n d e {l + V d ) 



(29) 



and is independent from y. 

Consider the case of underdense plasma n& > no, when all plasma electrons behind the 
bunch are "blow out" n d = 0,i.e. tp = —n ||14j|. Following |Tjj] assume that E d = i.e. 
if' Q = according to (|26|j27|) .Then constant h ( pl~l) is equal to 



TIq o 1 S -^S 7 

& = y + 3<^o = g n o = K 
i.e. the solution, corresponding to the "blow out" regime, lies on separatrix. 



The constant h (|21|),also can be expressed through the roots q = c^no of the equation 



1 aia 2 «3 3 
= -cic 2 c 3 = n 



(30) 



where —3/2 < ai < —1,-1 < a 2 < 0, < 03 < 1/2 for A; > 0. For the separatrix 
cti = «2 = —1, «3 = 1/2 



fl = lie 



^0 
6 



(31) 



For the values h > h s and h < as it is evident, the solutions for (p 2 have an infinite 
values. When ip Q < c 2 the solution became unphysical,even for < h < h s 

5 FINITE NONLINEAR SOLUTIONS 



First consider the case when c 2 < ipo, ¥2 < C3. From general solution ([2C|), using known 



expressions for the elliptic integrals and elliptic functions [ITB].[|T6[,we have 



(f2{z) = c 3 - (c 3 - c 2 )sn 2 z 1 , 



(32) 



where 
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(33) 



7o = arcsm 



I C3 - ipo 
C3-C2' 



q 



'C3-C2 a 3 -a 2 



' c 3 - ci v « 3 - ai 
and -F(7o, 9) is the elliptic integral of the first kind,sn2;i- elliptic function. 
Using (D,(H),(|3) and ©-(H) it is possible to obtain 

d 

E y = -— = -y<p 2 (z) = -yn [a 3 - (a 3 - ai)sn 2 zi] 



(34) 
(35) 



(36) 



E 7 



-n (a 3 - ati) 



1/2 



1 (y 2 + \ ( ^ 



-n (a 3 - ai) 



1/2 



sn2zi 



a 2 — ai * 4 

1 — I sn Z\ 

a 3 — «i 



The length A n of the nonlinear wave is given by 



dcp 2 



v Jc 2 ^2[h - F( V2 )] 1/2 



u \2n J (q, 3 + 
xnl (a 3 + \a 2 \ 



1/2 



x 



2 \ a 3 + a]. 



In the linear case 



(37) 



(38) 



h — > 0, «2 — ► 0a 3 — >■ 0, «i — > — 3/2 

and A n -> A p = 

From (|3^) it follows, that is by the order of magnitude equal to E y ~ yno,or in the 
ordinary units 

E y ~ = 4nen p y, 

e v 

which coincides with the field inside the flat bunch uniformly charged with the plasma 
electron density n p . 
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From (R7|) the longitudinal component of the electric field is by order of magnitude 
equal to E z ~ (H^ + aj n^ 2 , or in ordinary units 

e^oV 2 / fo \2 / 

where a = ^f, E z (y, z = —d) = E z0 + E z2 y 2 in ordinary units, so at some conditions it 
can be larger than E y . 

In the case h = h s the changes due to nonlinarity are more drastic. 

As we have seen,this case corresponds to the conditions 

n b > no, n e (z = -d,y = 0) = n d = 0, 
and E zQ = 0, which resembles the blow-out regime in underdense plasma (|6]). In this case 



ci = c 2 



— ^0; c 3 — * c s 



-no, 



1 

3< 



h - F((p 2 ) -> ^(v? 2 + n f I — - (p 2 



no 
2 



From ( p0|) with the minus sign in the front of integral (tp 2 decreases,when z increases 



from — oo up to —d), using |T5[ it follows 



V?2 = -no + 2 



(39) 



where 



ip = 9 — y/n~o~(z + d),9 = In 



-^ 2 ) 1/2 +(§no) 



(f -^) 1/2 -(|n ) 

In the ( f40"|) it is necessary to use <^o = — + e, e > Oand pass to limit e — > at fixed ^. 
Then at ^ = 0, z — > — oo,and at ip — > — oo, z — > — d. 

The electric field components in the considered case of the blow-out regime are 



(40) 



E y = -yn (3th 2 ij/2 - 1) 



E z = 3n 



:; 2 /V Ko^ 

1 2 n , 



[l-th 2 i)/2)thi)/2 



(41) 
(42) 
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The maximum value of the transversal component of the field by the order of magni- 
tude is E™ ax ~ yn . 

Longitudinal component E z = 0,when z — > — d, (ip — > — oo) and when z — > — oo, {ip — > 
0). Its maximum value is at ipQ = — 1, 32 and by the order of magnitude is equal 

V 2 n °) 

Hence the maximum values of the field coincides with that in case when < h < h s .The 
difference is in the form of the wave:when h < h s the wake wave (|36| - |3~TD is cnoidal and 



when h = h s (blow out regime) the wave (f41| - [42j ) is solitary one. 



6 APPROXIMATE SOLUTION OBTAINED BY MULTIPLE SCALES METHOD 

In order to obtain the solution of the basic eq. (|T5|) valid for larger values of y it is neces- 
sary to solve (|15|) by some other methods.lt seems that the multiple scales approximate 
method (development of the derivative is suitable for this propose. 

The small parameter in question for considered case is e = ( ei f m ^ x j <^ 1 ; (y/ <^l,ip 



here is in ordinary units). According to [0 introduce the different scales variables in z 

z = z,z n = e n (z + d),n= 1, 2, 3 . . . 
and perform the following developments 

^(y,z) = (p{y,z 0} z 1 ,z 2 ,. . .) = (43) 

= e<Pi(y, z, z u Z2, ■■■) + e 2 (p 2 (y, z, z u z 2 , . . .) + e 3 cp 3 (y, z, z 1} z 2 , . . .) 

d d d 2 d 
d-z = dz- +e dz- 1 +e d^ 2 + --- 

dz 2 dzl dz dz 1 \ dz dzi dz\ j 

(For dimensionless function and arguments it is necessary to put e = 1 in the final results.) 

Substitution of the developments ([43]) in (|i~5l) gives the following set of the equations 
for subsequent approximations: 



12 



^ + rwi = (44) 
9 2 <P2 d 2 cp 1 #Vi 

-dzf +n ^ = - Vi w 2 d^; (45) 

(9^0 no<;/93 ^2 ^ 2 ^2 dz dzi \ dz dz 2 dzf J < ^ >1 

The general solution of eq. (H4T) is 



<Pi = a( Vl , z 1} z 2 )e~ W° + a*( yi , z 1} * 2 )e + W° (47) 



and y-dependence of the solution (47) comes from boundary conditions at zq = —di, Z\ 
0,z 2 = O.If 

dipi(y,z = -d,Z! = 0,z 2 = 0) 
E z = ~ = 



E v = Wy f(y) 



then 



Rea(y, z = —d, z\ = 0, z 2 = 0) = (48) 

= ~\ I (f(y)cosn d+ —=g(y) sin ^/n^d] dy 
2Jo\ v /rio / 

Ima(y, z = —d, Z\ = 0, z 2 = 0) = 

= 7,1 ( f(y)smy/n^d+ —=g(y) cos y/n^d] dy 
2 Jo \ v m / 

The second term is right hand side of the eq. (fi~5|) is secular due to the solution ([47l);it 
can be eliminated if ipi is independent on Z\ i.e. 

It means that <p 2 is also independent on z\ and the general solution of the eq. 
(H3), taking into account (|49D,is 



<p 2 = B(y, z 2 ) + A(y, z 2 )e~ 2l ^ zo + A*(y, z 2 )e 2i ^ Z0 + b(y, zje'^* + b*(y, z 2 )e* 

(50) 
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where 

B = 1 (a d2a * | a* d2a ] A- —a— (51) 
no \ dy 2 dy 2 J ' 3n dy 2 



Function b(y,z 2 ) entering in the solution of the homogenious part of the eq.(K5|) can 
be found from the boundary conditions 

¥>2{y,Zo = ~d,z 2 = 0) = 

dip 2 (y,z = -d,z 2 = 0) 
dz 

and has the following value 







b = 1 (A*e- 3i ^> d - 3Ae l ^ d + B) (52) 

The 2 2 -dependense of the function a(y, z 2 ) (and subsequently the z 2 -dependense of 
A, B, b) comes out from the consideration of the eq. ([5]) for the third approximation. Eq. 
(f46[),due to the independense of <p±, (p 2 on z 2 (^9|) ,has the form: 

dVa j ^_ 2 gVi y^Va y/ 2y?1 (53) 
dz dz 2 dy 2 dy 2 

Using solutions (^7| , |50|) for ipi and ip 2 it is evident that right hand side of the eq. ([53]) 
has the secular terms, proportional to e ±l ^™° z ° .The conditions for their elimination are 

^ l_&af,d*a 0V\ 1 d 2 a d 2 a* 
dz 2 no dy 2 \ dy 2 dy 2 J 3n dy 2 dy 2 

—a— ( a *— + a—) —a*— (a—] 
n o dy 2 \ dy 2 dy 2 J 3n dy 2 \ dy 2 J 

and subsequent conjugate expression. 

The eq. ( |54"D determines the dependense of a(y,z 2 ) from z 2 .Eq. ( |54"D is complicated 
enough, it is a system of the first order differential equations for Rea(y, z 2 ) and Ima(y, z 2 ). 
The y-dependense of a(y, z 2 = 0) is given by boundary conditions. 

The eq. (|54|) simplifies under the assumption that a(y, z 2 ) is 

a(y,z 2 )=Y(y)Z(z 2 ) (55) 
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where Z{zi = 0) = l,and Y(y) is known from the boundary conditions (^8]) .Under (|55|) 
eq. (Q) takes the form 

dz 

2i^— = ^(y)\Z\ 2 Z (56) 

where if>(y) is 



m = -^—Y"Y*" * +^~^{Y" 2 - Y'Y'" - \yY"") 
3n 3n Y 2 

2 1 1 

_l _(Y Y f " + Y'Y* + — Y*Y"" + —YY* ) 

n 2 2 



(57) 



The function ip(y) is slowly varied on y,when \y\ < b and is zero,when \y\ > b.lt 
is reasonable to average the eq. fl56|) on y. Strikly speaking the need of such kind of 
procedure indicates that assumption ( 1551 ) is not in full appropriate, but practically can 
work for slowly varying ip{y). 

After averaging,eq. ( |56"D and its conjugate one give the following system for x = ReZ 
and y = imZ 



where 



(11* — — 

2J^— = Imi\)\Z?x + Reip\Z\ 2 y 
dz 2 

— 2^/nQ— — = Reijj\Z\ 2 x — Imif)\Z\ 2 y 

aZ2 



1 r b 

■0 = r / 

o JO 



(58) 
(59) 



From (|58| - |59D follows 



d|Z| 2 2Irmp lr/l , 



The eq. ( p0|) has the solution 



1 



c — 



z\\ 



> o, 



(60) 



(61) 



with the arbitrary constant c > ^=.In the simplest case, when Irmp = 0, \Z\ 2 = c 1 and 
the system (KSJ- |5D|) has the solution 



j; = c 1//2 cos 



y = — c *' 2 sin 



2 v /r^c 
' Reip 



Z2 + 0O 
Z2 + 0O 



(62) 



15 



The constants c 1 ^ 2 = 1,6q = 0,due to the boundary condition Z(z 2 = 0) = l.When 
Imip 7^ 0,x = \Z\ cos9,y = — \Z\ sin9, where 6 is the solution of equation 

d9 \Z\ 2 n 7 

Reip (63) 



dz 2 

at boundary condition Z(z 2 = 0) = 1, (6 = 0). 



» = \** H1 _*& H) (64) 

4 Imip Jtiq 



and a(y,Z2) for the first approximation fl4"T|),is a(y,Z2) = Y(y)Z(z 2 ) = Y{y)\Z\e %e , where 
\Z\ is given by (|6"ID,fl by (|64"D and Y(y) can be found from boundary condition (|48|). 



When Imip — > solution (|B1]) , (|64"|) turns to the solution (62). 



The solution fl6"4]) is valid for \z 2 \ = £ 2 \z'\ = ( e ^2 X ) k'l < 1 i>e. for = /c p |z| < 
f m ° 2 ~) .In the considered domain z 2 < 0,(z < —d) and the solutions (IBTl) , (|64T) have 
different behaviour, when Imip > and 7m?/> < O.When Imip > from (|6T|) it is seen 
that |Z| 2 (and consequently the amplitude a(y,z 2 )) decreases when \z 2 \ increases; when 
Imip < 0, \Z\ 2 and amplitude a(y,z 2 ) increases,when \z 2 \ = e 2 \z\ increases up to allowed 
value: 



mc \ / mc 



rain 



n _ I 



^^Pmax ) y&^Pmax J 2 / ITllb I 



An outlined approximate procedure, based on multiple scales method for solving basic 
nonlinear equation (|15|) for potential wake waves, shows that in the lowest (first, second 



and third) approximations solution of eq. (|T5|) is represented by set of harmonics ~ 
e iy/nZmz( m = o, 1, 2, 3 . . .) with amplitudes (a, B, A, b eqs.(|3|)2|), (0),(§J|3) slowly in- 
creasing or decreasing with |z|.The z-dependense of amplitudes is rather combersome.but 
it is clear,that some singularity can appear at certain boundary conditions when Imip < 
and I z 2 1 ~ TT^V < ( m ° 2 I , which is inside of the region of applicability of the adopted 
procedure. 

May be it is an indication of some new quality of the considered nonlinear poten- 
tial wake wave and, if it is so, it needs an additional consideration. In any case the 
oulined consideration indicates that boundary conditions at the rear end of the driving 
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bunch, which can be changed by appropriate choise of the bunch transverse and longitu- 
dinal charge distributions, can essentially effect the kind of the nonlinear wave amplitude 
dependense on the longitudinal coordinate z,even in the first approximation. 

It seems,that subsequent experimental investigation of the dependense of the wake 
wave amplitude on the driving bunch transverse and longitudinal charge distributions 
could be useful. 
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